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Abstract 

We discuss a new cluster representation for the internal energy and 
the specific heat of the d-dimensional Ising model , obtained by study- 
ing the percolation mapping of an Ising model with an arbitrary set of 
antiferromagnetic links. Such a representation relates the thermal op- 
erators to the topological properties of the Fortuin-Kasteleyn clusters 
of Ising percolation and is a powerful tool to get new exact relations 
on the topological structure of FK clusters of the Ising model defined 
on an arbitrary graph. 
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1 Introduction 



It is well known that the Ising model can be mapped into a percolation 
problem IQ, |2| . In its original formulation this mapping is based on the iden- 
tification of the mean value of the magnetization of the Ising model with the 
size of a (suitably defined) percolating cluster. Such identification is highly 
non-trivial and allows for a geometric characterization of several statistical 
observables near the critical point, a fact which has greatly improved our 
understanding of the Ising model and inspired new powerful algorithms ^] 
to simulate the model. This mapping can be extended also to the thermal 
sector of the model and cluster representations of, for instance, the internal 
energy and the specific heat can be easily constructed. However these rep- 
resentations, in contrast to the case of magnetization, are rather trivial and 
do not add any new geometrical or topological information to the model. 

In this paper we shall discuss a new cluster representation for the ther- 
mal sector of the Ising model. To this end we shall need first to extend 
the percolation mapping to the case of the Ising model frustrated by some 
antiferromagnetic links ||5| and then we use this mapping to extract a useful 
piece of information on the structure of the Fortuin Kasteleyn (FK) clusters 
of the unfrustrated model. 

The main feature of our new representation is that it relates the thermal 
operators of the model (and in particular the internal energy and the specific 
heat) to the topological properties of the clusters of a typical configuration 
of Ising percolation. As a consequence one can write down exact relations 
on the topological structure of the FK cluster. 

The simplest of these new relations deals with the concept of non-cutting 
or black bond. A bond of a FK cluster is said to be black if its cancellation 
does not split the cluster in two disjoint parts. We shall prove (see sect.^ 
that in the Ising model defined on an arbitrary d-dimensional lattice with 
A'^ links at the coupling /3, the mean number {Nb) of black bonds is related 
to the internal energy E by 

Combining this with the well known expression (see eq.(^) below) 

E = {Ng)-^-N , (2) 
smhp 

where (Ng) is the mean number of bonds of FK graphs associated the con- 
figurations of Ising percolation, yields a simple, exact relation between (Nb) 
and {Ng)- Similar relations can be found for the specific heat . 
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Actually these new relations are true for the Ising model defined on any 
arbitrary graph, however in extracting our results we shall always assume for 
sake of simplicity that the model be defined on a regular lattice in d = 2, 3, . . . 
dimensions with periodic boundary conditions. 

In the process of constructing our new representation we shall discuss 
two issues which are rather interesting in themselves. 

First, we shall construct a scheme to classify the bonds of a cluster on the 
basis of their topological properties (see sect. |3| below). The standard clas- 
sification Q, which splits the bonds in three classes (which are commonly 
denoted with the three colors red, blue and green) refers to the bonds of the 
(infinite) percolating cluster. Our new representation of thermal operators 
suggests a slightly different bond partition, based on the connection prop- 
erties of all the clusters. Within our classification scheme the topological 
properties of the bonds are unambiguously identified by the class to which 
they belong. 



Second, we shall show in sect. 4.1, using duality, that in the Ising model 



all the correlators of an even number of spins (and hence invariant under the 
Z2 symmetry of the model) can be expressed as ratios of partition functions 
with a suitable set of antiferromagnetic links. 

This paper is organized as follows. The first part of the paper is com- 
posed by an introductory section on the Ising model and the percolation 
map (sect. |2|) which will allow us to fix notations and to make the paper as 
self-contained as possible. In sect. ^ we discuss the topological properties 
of the clusters. In sect. Q we deal with the Ising model in presence of an- 
tiferromagnetic links and duality transformation and sect. |5| is devoted to 
the extension of the percolation mapping to the frustrated Ising model. In 
these last two sections we have collected most of our new results which are 
then used in sect. ^ to construct new cluster representations for the inter- 
nal energy and specific heat. Finally sect. |^ is devoted to some concluding 
remarks. 

2 Ising Model 

The d dimensional Ising model is defined by the Hamiltonian: 

H{J,h') = -J 'nSm + h'J2sn, (3) 

{n,m) 

where the field variable Sn takes the values —1 and -|-1; n labels the sites 
of the lattice (denoted with A in the following) which we assume to be a 
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d dimensional simple (liyper)cubic lattice of size L with periodic boundary 
conditions, however a large part of our considerations are valid for a Ising 
model define on an arbitrary graph. The notation (n, m) indicates that 
the sum is taken on nearest neighbour sites only. The partition function is 
defined as usual by 

Z= e-^^W'^') (4) 

Sn = ±l 

where f3 = Plugging eq. @ in the definition of Z and assuming the usual 
conventions: J = 1 and h = (3h' we obtain: 

Z{f5,h)= ^ + . (5) 

Sn=±l 

For /i = and d > 2 the phase diagram of the model is composed by two 
phases separated by a second order phase transition located at Pc ^ In 
the high temperature phase the Z2 symmetry of the model is preserved and 
the magnetization is zero, in the low temperature phase the Z2 symmetry is 
spontaneously broken and the magnetization becomes different from zero. 



2.1 Mapping to a percolation model 



In this section we shall discuss the mapping in the case h = 0. The extension 
to a non zero magnetic field is straightforward and can be found for instance 
in§. 

The Ising partition function of eq.(0) with h = can be rewritten as: 



ZiP) 



=±1 {ij) 



(6) 



where N = dL'^ is the number of links in the lattice and the 6 function takes 
the value S{si, Sj) = 1 when the two arguments coincide and zero otherwise. 
Expanding the products in eq.® one finds 



Z{(3) 



E E 

G Si=±l 



n p^'^^i 



(l-p) 



N-N(G) 



(7) 



where p = 1 — e^'^^, G denotes an arbitrary subgraph of the lattice and 
N[G) is the number of links of G. In general G will be composed by several 
connected components ( FK clusters in the following). Let us call G{G) the 
number of clusters in the graph G. Notice that among the clusters one has 
to consider also those with one site only. 
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Summing on the spin configurations in eq.(|^) we end up with 
= e^^^ _p)iV-7V(G) 2C(G) 

G 

which can be interpreted as the partition function of a percolation model 
with bond probability p and with a weight 2 for each independent cluster. 
In this framework the magnetization transition of the Ising model becomes 
a percolation transition, located at the percolation threshold = 1 — e~^^= . 

For p > pc an infinite, percolating cluster exists. The density of sites 
belonging to this percolating cluster, which is zero below p^, can be used 
as order parameter for the percolation transition. It can be shown that it 
exactly coincides with the magnetization density of the original Ising model. 

The internal energy in the percolation framework can be constructed by 
taking the logarithmic derivative of eq.(^) with respect to /?. The result is 

p 

where the mean value is taken with respect to the measure of eq.(Q). It is 
interesting to compare this result with the standard definition of internal 
energy: 

E = ^logZ{(3, h = 0) = {J2 SnSm) (10) 

(n,m) 

If we denote with the number of links which join spins with the same 
sign in a given configuration and = N — the number of those which 
join spins with opposite sign then eq.([To|) can be rewritten as: 

E = {N+)-{N^)=2{N+)-N . (11) 

By comparing eq.s (^) and ( |Tl[ ) we see that 

p{N+) = {N{G)) , (12) 

which gives the most intuitive way to define Ising percolation (and inspired 
Swendsen and Wang in their proposal [^): given a generic configuration of 
the Ising model, delete all the bonds which join spins with opposite sign 
and, on the remaining graph, construct a standard percolation process i.e. 
switch on the bonds at random with probability p. The resulting graph G 
will be a typical configuration of an Ising percolation model. 
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3 Cluster structure 



As we have seen, an interesting feature of the percolation mapping is that it 
allows a geometric characterization of various thermodynamic quantities. In 
order to better understand this geometric setting it is convenient to study in 
a more precise and detailed way the cluster structure in a typical percolation 
configuration. 

An important step in this direction was made by Stanley in 1977 |^ who 
noted that at p = pc, in a generic percolating cluster we can distinguish 
three different sets of bonds. By associating an electric unit resistance to 
each bond, and applying a voltage between the ends of the cluster, one can 
select the "dangling bonds" (also called green bonds) which are those which 
do not carry current. The remaining bonds form the backbone; in this set 
one can then select the singly connected bonds (red bonds) , which carry the 
whole current and have the property that if one is cut then the cluster breaks 
in two parts. The remaining bonds are multiply connected and are usually 
denoted as blue bonds. Starting from this coarse grained classification, one 
can then look at more subtle structures, selecting for instance pairs of double 
connected bonds, triples... Equivalently one can look to the subsets of bonds 
which carry exactly half of the whole current, one third, and so on. 



la) Stanley classification of the bonds of 
a percolating cluster. The thin lines de- 
note the dangling bonds, the thick lines 
the backbone and the dashed lines the red 
bonds. Free b.c. are understood 



lb) The same graph of fig. la in the topo- 
logical classification described in the text. 
Thin lines are the grey bonds, thick lines 
the black bonds and the dashed lines the 
red bonds. Periodic b.c. are understood. 
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With respect to this standard classification, which assumes free boundary 
conditions, in the present work we prefer to deal with periodic boundary 
conditions. This is because, as we shall see, the new representation of the 
thermal operators does not depend explicitly on the percolating property of 
the configurations, but rather on their winding and connection properties. 
We are then forced to introduce a different partition of the bonds which 
slightly differs from that of Stanley and applies not only to the (infinite) 
cluster wrapped along one or more directions, but also to any FK cluster 
contributing to the partition function. 

We suggest to the reader to follow our definitions looking at fig.l, where 
the bonds of a percolating cluster are split into three sets according to the 
Stanley classification (fig. la) and ours (fig. 2b). 

G denotes the generic graph which appears in the sum of eq. (^). It is 
composed by disconnected FK clusters. We call Gi the set of those bonds 
which, when deleted, disconnect the cluster to which they belong into two 
disjoint clusters. Remember that we consider as a cluster also that composed 
by only one site and no bonds. We call the bonds of Gi grey or cutting bonds. 

Let us call Go the complement of Gi in G, i.e. G = Go U Gi. Go 
contains two kinds of bonds: those which, if deleted, do not change the 
wrapping properties of G are called hlack bonds; those which transform a 
wrapped graph into an unwrapped one are called red bonds. 

It is clear that our grey bonds form a subset of the dangling bonds of 
the Stanley's classification, but there are dangling bonds that are black in 
our classification, as comparison of fig.s la) and lb) shows. 

Using again an electric circuit analogy, one can say that our black and 
red bonds are those which carry current in presence of a variable magnetic 
field pointing in a generic direction ^ 

It is known that the winding properties of the FK clusters are related 
to the interface tension in the low-temperature phase |1C]. In this paper we 
shall see that the black bonds are directly involved in a new representation 
of the internal energy. 

Our classification can be easily iterated by identifying other topological 
subsets of Go associated to other thermal observables. We need in particular 
to define another set, denoted with H2, that will play a major role in dis- 
cussing the new representation of the specific heat. H2 is the set of "cutting 
pairs" of Go, i.e. the set of those pairs of black bonds which, if simultane- 
ously deleted, disconnect the cluster to which they belong into two separate 



^We thank Lev Shchur for this observation. 
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Let us see a specific construction of the set H2, whicli is interesting 
both because it sheds hght on the structure of H2 and because it can be 
straightforwardly iterated to construct "cutting triplets" and from them 
possibly further higher combinations which contribute to higher derivatives 
of the free energy (in this paper however we shall not study these more 
complicated cases): 

• Choose a link I € Go and construct the graph Gq{1) obtained from Go 
by deleting the link I. 

• select in Go(0 the cutting bonds (as we did above when we constructed 
Gi starting from G), which form a set Gi(/). 

• repeat the operation for all / € Go- The set of pairs {(Z,/i),V/i € 
Gi{l),\/l G Go} is exactly twice H2 (due to the fact that each pair 
appears two times). This means that (apart from a factor two) H2 
coincides with the collection of graphs Gi{l) . 



4 Ising model with a set of antiferromagnetic links 

A natural extension of the Ising model discussed in sect. |2| is the one in 
which we allow the coupling J to change sign from link to link. Those 
links for which J = — 1 are said to be antiferromagnetic or flipped links. In 
analogy with the definitions of sect. |^ (and setting for simplicity h = 0) we 
have now 

H{J) = - ^ J{n,m)SnSm , (13) 
{n,m) 

from which we may construct the partition function: 

Z{13,J)= Y: e-^'^'^ (14) 

Sn=±l 

We denote the partition function without antiferromagnetic links and with 
periodic boundary conditions by Z{(3). Denoting by F the set of flipped 
links we can write 

= ( n e-^^-^^) , (15) 

^^^^ {^J)eF 



^Notice that, since by definition Go does not contain cutting bonds, at least two bonds 
are needed to disconnect a cluster. 
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where the mean values is taken with respect to the standard, unfrustrated 
Ising Hamiltonian. It is almost evident that suitable linear combinations 
of these expectation values reproduce the whole even sector of the theory. 
This can be simply proved using the Kramers -Wannier duality: This trans- 
formation gives a one-to-one map between these expectation values and the 
set of all correlation functions among the physical observables of the even 
sector of dual theory, as described in the next subsection. 

4.1 Duality and frustrations 

The Ising model defined in section ^ is characterized by the two follow- 
ing features: the dynamical variables (the spins) live on the sites and the 
interaction is defined on the links of the lattice A. Sites and links are respec- 
tively and 1 dimensional simplexes of the lattice. We can easily generalize 
the Ising model by looking to different geometrical realizations of dynami- 
cal variables and interaction Hamiltonians. For instance, in the Ising gauge 
model we take the spins si to live on links / and the interaction to be defined 
on plaquettes sq = IlieD • could also choose the dynamical variables 
on the plaquettes and the Hamiltonian to live on cubes, thus defining a 
Kalb-Ramond type theory. 

It is well known that the Kramers -Wannier duality transformation can 
be generalized to Ising models in any dimension d. This transformation is 
not a symmetry of the model: it maps one description of the dynamical 
system to another description of the same system. This duality transforms 
a given lattice A (in our case the (hyper)cubic lattice on which the Ising 
model is defined), into a new lattice (the dual lattice A) in which each k 
dimensional simplex is mapped to a, d — k one. In the case of (hyper)cubic 
lattices the dual lattice turns out to be again (hyper)cubic. 

Under this mapping the Ising model on A is transformed into a new 
model whose Hamiltonian lives on a d — 1 simplex of A. Thus we have the 
following correspondences: 

if d = 2 Ising model <^=^ Ising model 

if d = 3 Ising model <J=^ Z2 gauge model 

if d = 4 Ising model<^=^ Z2 Kalb-Ramond model. 

It can be shown that in the thermodynamic limit the free energy density 
of the original Ising model coincides (apart from an additive function of /3 
which can be evaluated exactly) with that of the dual model evaluated at 
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the dual coupling (3 defined as 

^ = _i \og[thm . (16) 

Using a duality transformation it is possible to build up a one-to-one map- 
ping between a given pattern J = {J{ij)} of antiferromagnetic links and 
the correlators of the physical observables of the dual description. Take for 
instance the 3-D case and denote by F the set of flipped links. Associate 
to each of these links the corresponding plaquette of the dual lattice and 
denote by F this set. Duality implies that ^ 

( n e-2/5^'^^-) = (n 5D)gauge , (17) 

where the left-hand side coincides with the ratio of the partition functions 
defined in eg. (p^) , while the right-hand side is the expectation value of the 
product of plaquette variables sq with a Boltzmann factor e~^'^saugo -with 

•Sgauge = — ^ SQ ; SQ = Sl^Sl^Sl^Sl^, (18) 

□ 

where s/ G {1,-1} are the variables located in the links of A. 

Conversely, given any product of Wilson loops W{C\)W{C2) . . ., we can 
replace it (in many ways) with an equivalent set of elementary plaquettes, 
with the only constraint that the boundary of this set should coincide with 
the set of loops {d,C2, ...}, then we can apply again eq.(p!7|) and this 
completes the proof of the above-mentioned one-to-one correspondence in 
the 2>D case. The extension to other dimensions is straightforward. 

We want now to use this correspondence to show that for the even sector 
of the spin Ising model (namely for the operators obtained as products of 
an even number of spins and hence invariant under the Z2 symmetry of the 
model) the set of ratios (|T5|) corresponding to all the possible choices of an- 
tiferromagnetic bonds forms a complete set in the sense that any correlator 
of the spin Ising model can be expressed in terms of these ratios. To be 
definite we put forward this argument again for the 2>D case, its general- 
ization to other dimensions being straightforward. The set of products of 
Wilson loops forms a complete set of observables of the dual description; the 

^It is important to stress that eq.(]l^) holds only in the thermodynamic limit, in which 
the boundary conditions can be neglected. In fact in general the duality transformation 
does not preserve specific choices of the boundary conditions. In particular, periodic b.c. 
are mapped by duality into a mixture of periodic and antiperiodic boundary conditions. 
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simplest, non-trivial observable is the plaquette. The spin-spin correlators 
on the original lattice A can be expressed in terms of the gauge model by 
flipping the set of plaquettes of A intersected by a path 7 C A connecting 
the site x to the site y: 

(S.^.) = (ne"''^'°)sauge (19) 

Using the obvious identity e~'^^'^^ = cosh 2/3 — sinh2/3sQ the right-hand 
side can be written as a combination of expectation values of products of 
Wilson loops. Then using eq.(|l7|) we can evaluate this correlator directly in 
the Ising model in terms of flipped partition functions. Such a construction 
can be repeated for any other spin correlator and easily generalized to any 
dimension d. In conclusion, we can state that the whole set of spin correla- 
tors of the Ising model in any dimension can be encoded in a suitable set of 
ratios of flipped partition functions. 



5 Percolation mapping of the frustrated model 



The percolation mapping discussed in sect. |2.l| can be extended also to 
the frustrated Ising model |^] ^ The main difference with respect to the 
unfrustrated case is that now also the topology of the cluster is important. 
Let us see this mapping in detail. 

We assume that the set of couplings Jij of eq.(^) is given by an arbitrary, 
fixed collection of signs ±1. 

As in the case without frustrations the Ising partition function of eq. ( p!^ ) 
can be rewritten as: 



Z(/3,J) = e^^ Y{[l-p + p5{J^i^)S,Sj) 



(20) 



where 6{J(^ij^SiSj) = (1 -|- J(^ij^SiSj)/2 is a projector on configurations with 



J(ij)SiSj = 1 . Expanding the products in eq. (P0[) we find: 



Z{P, J) = e^^ ^(1 - pf-^(G)p^(G) ^ 



G 



Si=±l 



(211 



^See also ||9| for some recent application of the percolation mapping in presence of 
frustrations to the study of disordered systems. 
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as in sect. |2|, is the total number of links of the lattice, the summation 
goes over all the subgraphs G of the lattice, N{G) is the number of links of 
G andp = l-e-2/5. 

Note that the whole dependence on the signs J^j^^ is contained in the 
6 projectors and that the product inside the square brackets is a projector 
which forces all the links of each connected component Gc of the graph to 
fulfill the constraint 

%)SiSj = 1 , {ij) G Gc . (22) 

We say that a cluster Gc is compatible with a given choice of couplings 
J^jj) if there is a configuration of its sites obeying such a constraint. In the 
standard case {J{ij) = 1 V (ij)) this condition implies simply that the sites 
of each cluster have the same sign. In presence of frustrations it is easy to 
verify that eq.(|2^) is a topological constraint, telling us that a connected 
graph Gc is compatible if and only if no loop of Gc includes an odd number 
of antiferromagnetic links. 

A crucial observation is that, owing to the connected nature of the clus- 
ter, the value of Sj of any site i of a compatible cluster fixes the values of all 
the other sites of the cluster, and that if {si , i £ Gc} is a solution of the 
constraint, also the opposite configuration {— Sj} is a solution. Thus any 
connected subgraph contributes to the partition function with a factor of 
2 if it is a compatible cluster, otherwise it gives a zero contribution. Thus 
summing on the spin configurations in eq.(|2l|) we end up with 



ZiP, J) = e^^ J2 ^XG')p^(^) (1 - 2^(^) . (23) 

G 

where zuj denotes the projector on compatible graphs, i.e. graphs made with 
compatible clusters, defined as follows 

1 if no loop of G contains an odd 
wj(G) = ^ number of antiferromagnetic links (24) 
otherwise . 

In the standard, unfrustrated case the sum over G is unconstrained and 



we obtain the result discussed in sect. 2.1. When there are frustrations, the 



set of compatible graphs is a proper subset of the all the possible subgraphs 
of the lattice. We can then write the following exact relation, which is the 



cluster version of the eq.(15): 

ZiP.J) 



ZW) ■ 
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where the expectation value is taken with respect the standard Hamiltonian. 

An interesting, particular example is the Ising model on a cubic lattice 
with periodic (p) boundary conditions in one coordinate direction, say z, in 
which all the links of a slice orthogonal to z are flipped to —1. Such a pattern 
of antiferromagnetic couplings is equivalent to choose antiperiodic b.c. (a) 
along the z direction. Denoting with Za and Zp the partition functions with 
these two different choices of b.c, we have 

^ = {^z) , (26) 

where Wz is a projector on the FK graphs which is 1 if there is no cluster 
with an odd winding number in the z direction and otherwise. Eg. ([2^) is 
the starting point of a new representation of the interface free energy first 
found in ref. which is particularly useful near the critical point. 

Coming back to the general formula (pSj), we would like to make a few 
comments: 

• Eq.(^) is valid not only for any regular lattice in any dimension, but 
also for the Ising model defined on an arbitrary graph and for any 
choice of antiferromagnetic couplings. 

• Any correlator of the even sector of the Ising model can be expressed 
in terms of the ratios (p5|). 

• The projector w j depends only on the subgraph Go C G. 

As a consequence, there is no information loss in the even sector if one delete 
all the grey bonds of any graph G, provided one use eq.(|^) to evaluate these 
observables. In other words, in the partition function (P) we can split the 
sum over all possible graphs X^g double sum 

ZiP) = e^PJ2 E p^(^)(l-p)^-^(^)2^(G)=^g-7Y(Go;/3)^ (27) 
Go Gi(Go) Go 

where G = Go U Gi and Gi(Go) denotes an arbitrary set of grey bonds 
compatible with a fixed set Go of black (and red) bonds. Tl defines a new 
Hamiltonian which depend only on the configurations of the black bonds. 



Eq. (25) tell us that, in spite of the sum over all the possible insertions of 
grey bonds, in the resulting Hamiltonian TC it is encoded exactly the same 
piece of information about the even sector than that of the original Ising 
Hamiltonian. Though this fact will not be exploited in the present paper, it 
clearly suggests the existence of an hidden huge symmetry of the theory. 
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6 Cluster description of thermal observables 



We are at this point in the position to study a new cluster representation 
of thermal observables, alternative to the one presented at the end of sect. 



2.1. Let us look first at the internal energy. 



6.1 Internal energy 

Let us consider an Ising system defined on an arbitrary lattice in d space 
dimensions with only one antiferromagnetic link in the position (ij). Then 
eq.(^) can be rewritten explicitly as 

^^^'/^y^^ = = cosh 2/3 - (siSj) smh2(3 . (28) 

Z{(j) 



On the other hand, using eq.(25) we get 



cotanh 2/3 -^M^ , (29) 

where ro^jj^ is the projector on the graphs compatible with the antiferro- 
magnetic link in the position (ij). Such graphs are of two types: 

• those in which the link {ij) does not appear; 

• those in which the cluster which contains the bond (ij) is split into 
two separate clusters when it is deleted. 

These two conditions tell us simply that (ij) does not belong to the subgraph 
Go of G. Summing over all the possible links and taking advantage of the 
translational invariance of the lattice we find: 

(E ^{nm)) = N - {NiGo)) . (30) 

(n,m) 



Combining this result with eq.s(29) and (|lO|) we get the sought-after repre^ 



sentation of the energy in terms of black bonds 

i? = ^^^^ + A^tanh/3, (31) 
smh(2/3) 

as anticipated in the introduction. This is the first main outcome of our 
analysis. 
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6.2 Specific heat 

Let us start from the following definition of the specific heat 

^ = ((-Yl SnSm){Y SkSi)) - SnSm)f (32) 

(n,m) {k,l) (n,m) 

call / = (nm) and g = (kl) the two links and separate the case f = g in 
the sum. We obtain 

C={N-^) + Yl ii'f's) - {sf){^a)] > (33) 

with Sf = SmSn Stnd Sg = SkSi. 

In order to get the new cluster representation of C it is sufficient to 
consider now a system with two antiferromagnetic links located in / and g 
with the associated projector TUfg. Then we have 

{wfg) = (e-2^("/+^9)) = ((cosh2/3-s/sinh2/?)(cosh2/3-Sgsinh2/3)) . (34) 

Using again eq.(|29|) we get 

EK^/-.>M..)M^E '"^t^r/;<--' . (35) 

f¥=9 i+9 ' 

The graphs with ■uofg = \ are of two types: 

• those in which both / and g do not belong to the Go subgraph; when 
summing over / and g these configurations give simply the contribu- 
tion E/g vofWg = {N- N{Go)){N - N{Go) - 1). 

• those in which both / and g belong to Go but are such that any loop 
going through / contains also g. Let us denote by Gjg the projector 
which selects the graphs with this special property. 

We may give a cluster representation of the specific heat if we are able to 
evaluate the sum Ylf^g Gfg- Reversing the order of summations, i.e. taking 
first the sum over all the pair of links and then the sum over G implied in 
the expectation value we obtain 

Y,{gfg) = 2{N{H2)) (36) 

In fact, from the definition ol Gfg, ii Gfg = 1 then cutting simultaneously / 
and g we disconnect the cluster to which they belong into two subclusters. 
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and conversely if two links belong to Go and are a cutting pair then they 
certainly fulfill the condition Gfg = 1 . The factor two is due to the fact that 
in the sum the pair /, g appears twice. 

Collecting together the various pieces we end up with 



which gives the new representation of the specific heat in terms of black 
bonds. 

7 Conclusions 

In this paper we have shown how to construct a new representation of ther- 
mal operators in Ising percolation in terms of a new set of bonds, called the 
black bonds, forming a subset of the standard FK clusters. Our approach is 
quite general, indeed we never needed to specify the lattice structure. Ac- 
tually our results are true also for Ising models defined on arbitrary graphs. 
Our main results are: 

• We proposed a new scheme to classify the bonds of a cluster on the 
basis of their topological properties. We have discussed the relations 
with the standard classification scheme. 

• We pointed out that in the Ising model all the correlators of an even 
number of spins can be expressed as ratios of partition functions with 
flipped links. In particular we have applied this result to the internal 
energy and the specific heat. 

• We found that in the configurations of black bonds contributing to the 
partition function is encoded the whole information on the even sector 
of the theory; in particular the internal energy can be expressed in 
terms of the mean number of black bonds and that the specific heat is 
related to their variance and to the special subset H2 defined in sect. 



• finally we observed that our new representation of observables of the 
even sector suggests the existence of an hidden huge symmetry of the 
Ising model. 

Our results can be straightforwardly extended to g-state Potts models 




sinh2(2/3) 



(37) 



i 



with generic values of q fl^ . 



15 



The most relevant application of the present analysis is that it may 
help to find new powerful algorithms to simulate the Ising model or to con- 
struct improved estimators for thermal observables. As a matter of fact 
for some special configurations of antiferromagnetic links these algorithms 
already exist. In the case of a whole flipped hyperplane leading to antiperi- 
odic boundary conditions this possibility was discussed for the first time by 
M.Hasenbusch in [|lO|, who provided a very powerful tool to evaluate the 
surface tension. Later, it was modified so as to evaluate in the dual version 
Wilson loops [0], correlators of Polyakov loops [|l2| and plaquette expecta- 
tion values ||l3[ in the gauge Ising model. We hope that the present analysis 
could help to further extend the range of these applications. 
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